Discrete Cosine Transforms on Quantum Computers 
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Abstract 

A classical computer does not allow to calculate a dis- 
crete cosine transform on N points in less than linear time. 
This trivial lower bound is no longer valid for a computer 
that takes advantage of quantum mechanical superposition, 
entanglement, and interference principles. In fact, we show 
that it is possible to realize the discrete cosine transforms 
and the discrete sine transforms of size N x N and types 
I, II, III, and IV with as little as 0{\og^N) operations on a 
quantum computer, whereas the known fast algorithms on a 
classical computer need 0{N log N) operations. 



1 Introduction 

Feynman proposed in 1982 a computational model that 
was based on the principles of quantum physics instead of 
classical physics. The model has been considered a mere 
curiosity until Peter Shor showed in 1994 that it is possible 
to factor integers in polynomial time on a quantum com- 
puter Thus, a moderate sized quantum computer is for 
instance able to break the RSA public key cryptosystem. 
Quantum computing is an exciting area of emerging signal 
processing applications. In fact, signal processing methods 
play a key role in Shor's integer factoring algorithm and in 
many other quantum algorithms. 

A quantum computer is based on the concept of a quan- 
tum bit, just as a classical computer is based on the notion 
of a bit. A single quantum bit represents the state of a two- 
level quantum system such as a polarized photon or a spin- 
1/2 system. Unlike a classical computer, adding another 
quantum bit to the memory of a quantum computer will not 
increase the dimensionality of the state space by one but will 
double it, allowing for linear combinations of 2" different 
base states in the case of n quantum bits. 

A program on a quantum computer is composed of a se- 
quence of elementary 'gates', which perform simple uni- 
tary transforms (as explained in Section 3). In fact, many 



algorithms in quantum computing rely on the fast Fourier 
transforms, the Walsh-Hadamard transforms, or other uni- 
tary transforms well-known in signal processing. 

The purpose of this paper is to derive (extremely) fast 
quantum algorithms for the discrete cosine and sine trans- 
forms. These algorithms can be implemented on a num- 
ber of quantum computing technologies based on Raman- 
coupled low-energy states of trapped ions [^, ^, nuclear 
spins in silicon [^, electron spins in quantum dots Q|, 
atomic cavity quantum electrodynamics [^J, or linear op- 
tics [01 . Coherent control of up to four quantum bits has 
been demonstrated and much progress is expected in 
the near future. 

2 Definitions 

Recall the definitions of the discrete cosine trans- 
forms: 
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where ki := 1 for i = 1, . . . , iV — 1 and fco l/\/2. The 
numbers ki ensure that the transforms are orthogonal. The 
discrete sine transforms S]^, Sj^, S]^^, and S'jy' are defined 
accordingly, see [§ for details. Notice that C^^ (resp. 5]^^) 
is the transpose of C]^ (resp. Sj^), hence it suffices to derive 
circuits for the type 11 transforms. In the following, we con- 
tent ourselves to iV = 2", which is justified by the machine 
model introduced below. 
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3 Quantum Gates 

A quantum computer consists of a system of n two-level 
quantum systems each of which represents a quantum bit 
(qubits). The computational state space is given by . 
Denote an orthonormal basis of by \x) where x is an 
rt-bit integer. The state of the quantum computer can be 
manipulated by quantum gates. Two types of operations are 
considered as elementary: the single qubit operations and 
the controlled NOT operations . 

The single qubit operations are given by local unitary op- 
erations of the form l2n-t ^U^l2t-i, with U E U{2). The 
controlled NOT gate operates on two qubits. It negates the 
target qubit if and only if the control qubit is 1 . All opera- 
tions on a quantum computer can be built up from these el- 
ementary operations, i.e., they form a universal set of gates. 
Therefore a basic (and nontrivial) task is to find efficient 
factorizations for explicitly given transformations into ele- 
mentary gates 

There is a graphical notation for quantum gates that has 
been introduced by Feynman. Each line denotes a qubit 
with the most significant bit on top. The circuits are read 
from left to right like musical scores. The figure on the left 
shows a single qubit operation U ® \2- The figure on the 
right shows a controlled NOT operation, that is, the unitary 
transform |00) ^ |00), |01) >-> |11), |10) ^ |10), |11) ^ 
|01). 
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We refer the reader to for more information about 
quantum gates and the graphical notation. 

4 DCTandDSTof Typel 

We derive the circuits for the discrete sine and cosine 
transforms of type I all at once. Indeed, the DSTj and DCTj 
can be recovered from the DFT by a base change Jll]] 
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and Fat = -^[exp{2Tri kl/N)]k,i=Q..N^i with = -1 
denotes the DFT of length N. Since efficient quantum cir- 
cuits for the DFT are known [[l|], it remains to find an effi- 
cient implementation of the base change matrix T^. 

Denote the basis vectors of C-^"^ by \bx), where 6 is a 
single bit and x is an n-bit number The two's complement 
of an n-bit unsigned integer x is denoted by x', that is, x' = 
2" — X. The action of T/v can be described by 



Tat |0O) = |00) , Tn\0x) = ^^\0x 
Tn\10) = \10 
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for all integers x in the range 1 < x < 2", where = —1. 
Ignoring the two's complement in T^, we can define an 
operator D by 

D|00) = |00), D\0x) = ^\0x) + ^\lx), 
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for all integers x in the range 1 < a; < 2". This operator is 
essentially block diagonal and easy to implement by a single 
qubit operation, followed by a correction. Indeed, define the 



1 



1 i 

1 -i 

implementation of the operator D. 



matrix Bhy B = —= 
V2 



then Figure |I] gives an 
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Figure 1. Circuits for the matrix D and the 
permutation tt. 



Define tt to be the permutation given by a two's comple- 
ment conditioned on the most significant bit tt jOx) = jOx) 
and TT |la;) = |lx') for all n-bit integers x. It is clear that 
T/v = TiD. The circuit for the permutation tt is shown in 
Figure |l]. Here P„ denotes the map \x) ^ mod 2") 
on n qubits, see [|l2| for an implementation. 

Theorem 1 The discrete cosine transform C\ and the dis- 
crete sine transform S\q can be realized with 0{\og^N) 
elementary quantum gates; the quantum circuit for these 
transforms is shown in Figure |] 

Proof Let N = 2". We note that 0(log^ N) quanfiim gates 
are sufficient to realize the DFT of length 2N, see |jl|]. The 
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Figure 2. Complete quantum circuit for the DCTi 



permutation tt can be implemented with at most 0(log^ N) 
elementary gates. At most O(logA^) quantum gates are 
needed to realize the operator D. This shows that the DCTj 
and the DSTj can be realized with 0(log^ N) elementary 
quantum gates. The preceding discussion shows that Fig- 
ure I realizes the DCTj and DSTj . □ 

5 DCT and DST of Type IV 

The trigonometric transforms of type IV are derived 
from the DFT by 

e"/4^i?^ . . ^ ^iv ^ (_^)5iv. (2) 
Here Rn denotes the matrix 
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with Lo — exp(27riMA^). Equation (||) is a consequence of 
Theorem 3.19 in [ [ll[ ] obtained by complex conjugation. 

Theorem 2 The discrete cosine transform Cj^ and the dis- 
crete sine transform S^^ can be realized with 0(log^ N) 
elementary quantum gates; the quantum circuit for these 
transforms is shown in Figure ^ 

Proof. It remains to show that there exists an efficient quan- 
tum circuit for the matrix Rj^i in equation (||). A factoriza- 
tion of Rn can be obtained as follows. Denote by x the 
one's complement of an n-bit integer x. We define a per- 
mutation matrix tti by tti |0a;) = |0a;) and tti \1x) = 
for all integers x in the range of < a; < 2". Denote by Di 
the diagonal matrix 



Di = diag(l, uj,..., cJ^, ...,uj\Ld). 

Then Rm can be factored as Rm = tti ■ Di ■ [B ® \n)- 



Note that i? ® Iat is a single qubit operation, and tti can 
be realized by controlled not operations. The implementa- 
tion of the diagonal matrix Di is more interesting. Note 
that 



Ai diag(l,LJ, 
A2 = diag(ZU^~ 
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where Lj — diag(l,a; ) and Kj — diag(a; ,1). 
Therefore, it is possible to write Di in the form Di = 
(C(g) Iat) • (Ai e A2) with C = diag(l,cJ). 

The complete quantum circuit for the DCTjy is shown 
in Figure |^. Note that the last three single qubit gates C, 
B^, and M ~ diag(e^'/"'^, e'^*/^^) can be combined into 
a single gate M B^C. □ 

6 DCT and DST of Type II 

The implementation of the trigonometric transforms of 
type II follows a similar pattern. Both transforms can be 
recovered from the DFT of length 2N after multiplication 
with certain sparse matrices, cf. Theorem 3.13 in 
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Figure 3. Complete quantum circuit for DCTiv 



and Lu = exp(27ri/4iV) with i'^ = —1. 

Theorem 3 The discrete cosine transform C]^ and the dis- 
crete sine transform S'JJ can be realized with 0(log^ A'^) 
elementary quantum gates; the quantum circuit for these 
transforms is shown in Figure Q 

Proof. We need to derive efficient quantum circuits for the 
matrices Vn and Un in equation The matrix Vn has a 
fairly simple decomposition in terms of quantum circuits. 

Lemma 4 Vn = i^i{H ®\im). 

Proof. It is clear that the Hadamard transform on the most 
significant bit 7f (g) Iat is - up to a permutation of rows - 
equivalent to Vat. The appropriate permutation of rows has 
been introduced in the previous section, namely tti |0a;) = 
I la;) and tti \1x) = \Vx) for all < x < 2". We can 
conclude that V^v = 7ri(-ff ® Ijv) as desired. □ 

The decomposition of is more elaborate. Notice that 
C/w|0O) = |0O) C/at 111) = (-1)110) 
Un\Qx) = ^|02;) + ^|1x') 
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for all integers x in the range 1 < x < 2" and all integers y 
in < ?/ < 2" — 1. Here and 1 denote the n-bit integers 
and 2" — 1 respectively. 

Define L»o by Do |10) = i |10) and L»o |a;) = \x) oth- 
erwise. We define a permutation 7r2 by 112 \^x) — |0x) 
and 712 \^x) — \l{x + 1 mod 2")) for all integers x in 
< X < 2". 

Lemma 5 C/at = _d| Tat _Dq 712. 

Proof Since d\ \Qx) = \0x) and d\ \lx) ^ uj""' \ lx), 
we obtain 

D\Tn\Qx)= ^\Qx) + ^\lx') 



D\Tn 
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We have Z)o7r2 \0x) = |0a;) and moreover Dott2 \1x) — 
|l(a; + 1 mod 2")) for all integers x in < a; < 2" - 1, 
andDo7r2|ll) = i |10)^ We note that (x+l mod 2")' = x, 
whence combining DiT^ with Dott2 shows the result. □ 
Recall that Tjq = ttD. It follows that 

Ul = 7r^\DoD')7r-'D,. 

The implementation of Di has been described in the sec- 
tion on the DCTiv, and the implementation of tt (and hence 
TT~^) is contained in the section on the DCTj. The imple- 
mentation of 772^ is also straightforward. It remains to find 
an implementation of DqD*^. We observe that 

DoZ?* |00) = |0O) , DoD'^ \0x) 
DqD^ 1 10) = -i 1 10) , DqD^ \lx) 
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This can be accomplished by a single qubit operation fol- 
lowed by a multiply conditioned gate, where the single 
qubit operation is given by i?* (g) Iat and the conditional 
gate acts via 
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The full circuit is shown in Figure ^. The statement about 
the complexity is clear □ 

7 Conclusions 

Signal processing methods have proved to be useful in 
virtually all known quantum algorithms. A basic problem 
in the design of quantum algorithms is to choose a well- 
adapted basis to recoup relevant information about the state 
of the system. The well-known decorrelation properties of 
DCTs may prove to be useful within this framework. We 
have shown that the DCT and DST of types I, II, III, and 
IV can be realized with a polylogarithmic number of ele- 
mentary operations on a quantum computer Compared to 
the classical realization of the DCT, this is a tremendous 
speed-up, making the DCT attractive in the design of other 
quantum algorithms. 
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Figure 4. Complete quantum circuit for DCTn 
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